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The low-frequency and shot noises in spin- valve CoFeB/MgO/CoFeB magnetic tunneling junctions 
were studied at low temperature. The measured 1/f noise around the magnetic hysteresis loops of 
the free layer indicates that the main origin of the 1/f noise is the magnetic fluctuation, which is 
discussed in terms of a fluctuation-dissipation relation. Random telegraph noise (RTN) is observed 
to be symmetrically enhanced in the hysteresis loop with regard to the two magnetic configurations. 
We found that this enhancement is caused by the fluctuation between two magnetic states in the 
free layer. Although the 1/f noise is almost independent of the magnetic configuration, the RTN is 
enhanced in the antiparallel configuration. These findings indicate the presence of spin-dependent 
activation of RTN. Shot noise reveals the spin-dependent coherent tunneling process via a crystalline 
MgO barrier. 

PACS numbers: 75.70.Cn, 73.50.Td, 73.40.Rw,72.25.Ba 



I. INTRODUCTION 

The magnetic tunneling junction (MTJ), which con- 
sists of a tunnel barrier sandwiched between two fer- 
romagnetic electrodes, is one of the central topics in 
spintronics. 1 MTJs exhibit tunneling magnetoresistance 
(TMR); their resistance depends on the relative magnetic 
configurations (parallel or antiparallel). Since the TMR 
effect was discovered by Julliere,— amorphous AI2O3 has 
been mainly used as a tunnel barrier 4^ However, in 2004, 
large magnetoresistance (MR) was obtained in MTJs 
with a crystalline MgO barrier-^ supported by theo- 
retical prediction^— These days, MgO-based MTJs are 
extensively studied from the viewpoints of fundamental 
physics and device applications. 

Although most MTJ studies have thus far been per- 
formed via conventional resistance measurements, noise 
measurements can serve to further clarify the intrinsic 
properties in MTJs. The noise results from the fluctu- 
ation of the current (thermal noise and shot noise) and 
of the resistance such as the 1/f noise and the random 
telegraph noise (RTN). Thermal noise and shot noise are 
due to the thermal agitation of electrons and the partition 
process of electrons, respectively, whereas the resistance 
fluctuation in MTJs is attributed to a nonmagnetic ori- 
gin (charge trap in the tunneling barrier) and a magnetic 
origin (magnetic fluctuations and domain wall motion in 
the free and/or fixed magnetic layers). 

Shot noise offers information on the interactions 
and/or quantum correlations of conducting elec- 
tronsj 10 ' 11 When the average current / is fed to a tun- 
nel junction, the current noise Si resulting from the 
shot noise can be expressed as Si — 2eIF (in the zero- 
temperature limit) with Fano factor F. It is well es- 
tablished that F = 1 in normal-insulator-normal junc- 
tions^ which means that the electron partition at the 
junction obeys a Poissonian process. In the MTJ case, 



when the tunnel barrier is composed of amorphous 
AI2O3, electron tunneling can be explained by the con- 
ventional Julliere's model. However, a coherent tunnel- 
ing via highly spin-polarized Ai Bloch state a 13 i 14 plays a 
central role in MTJs with crystalline MgO barriers. Al- 
though information obtained by conventional I-V mea- 
surements cannot directly address the coherence of elec- 
tron transport, shot noise offers further insight into the 
mechanism of electron transport J^r— In fact, we re- 
ported sub-Poissonian shot noise (F < 1), which is at- 
tributed to coherent tunneling, in a previous work^ Our 
experimental work is quantitatively reproduced by the re- 
cent theoretical work with first-principles calculations. 23 
Resistance fluctuations are also important, as MTJs 
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FIG. 1. (a) Schematics of measurement setup and device 
structure, (b) Typical MR curves of the present MTJ mea- 
sured at 5 K. (c) Typical differential resistance for P (dashed 
line) and AP (solid line) configurations as a function of bias 
voltage. 
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have found broad application, such as for magnetic field 
detectorsj2i magnetic random-access memory ; 25 ' 26 ran- 
dom number generators^ 6 -, and microwave oscillators.— 
For these applications, the signal-to-noise ratio is crit- 
ically important, where the 1/f noise and RTN limit 
device performance at low frequency. There have been 
many reports of 1// noise and RTN in MTJs with AI2O3- 
basedj2&~— MgO-based^ 5 . - — and other tunneling bar- 
riers i 22 ' 48 Nevertheless, little is known about the noise 
properties of MTJs with submicron-sized junctions with 
a thin tunneling barrier ; 46 ' 47 which are envisaged for 
memory and oscillator applications. 

In this paper, we report on noise properties, includ- 
ing shot noise, 1/f noise, and RTN in well-crystalline 
MgO-based MTJs with submicron-sized junctions with 
thin tunneling barriers. The noise measurement was car- 
ried out at low temperature with high experimental ac- 
curacy, with a focus on the noise properties around the 
magnetic hysteresis loops of the free layer. We investi- 
gated each noise source systematically as a function of 
magnetic field and bias voltage. The clear dependence 
of the 1/f noise on the applied magnetic field indicates 
that the main origin of the 1/f noise is magnetic fluc- 
tuation in the free layer, which is discussed in terms of 
a fluctuation-dissipation relation. Based on the bias de- 
pendence of RTN, we discuss the origin of RTN, which is 
different from that of the 1/f noise. The analysis of shot 
noise is also presented to further support our previous 
report^ - 

This paper is organized as follows: In Sec. IIA, 
we provide information on the sample fabrication and 
basic properties of our sample. Then the measure- 
ment system is described. In Sec. IIB, the anal- 
ysis method to extract the frequency-dependent and 
frequency-independent components is explained. Section 
III is devoted to the experimental results and discussion. 
First, we show the magnetic field and bias-current de- 
pendence of frequency-dependent noise in Sec. IIIA, and 
then the origin of the 1/f noise is discussed by using 
the fluctuation-dissipation relation in Sec. IIIB. We es- 
timate the fluctuating magnetic moment of the RTN in 
Sec. IIIC. The bias voltage dependence of the measured 
RTN is analyzed to obtain information on the excitation 
mechanism in Sec. HID, and, finally, we show the result 
of the white noise component in Sec. HIE. In Sec. IV, we 
conclude our study. 



II. EXPERIMENT 

A. Device and measurement 

Multilayer stacks of MTJs were deposited in a mag- 
netron sputtering system on a SiC>2 layer on a silicon 
substrate. The order of the layer structure from the 
substrate is as follows; buffer, PtMn(15), CoFe(2.5), 
Ru(0.85), CoFeB(3), MgO(1.05), CoFeB(2), and cap, 
where the top CoFeB layer serves as a free layer [see 
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FIG. 2. (a) Typical voltage noise spectral density of our MTJs 
in the AP configuration for different bias voltages: V = 12, 
22, and 55 mV at 5 K. The dashed line is a guide to the 
eye. The inset shows the same spectra between 160 and 200 
kHz (but note that the bottom axis is a linear scale), (b) 
Measured real-time voltage fluctuations (with vertical offset 
for clarity) for the same condition as (a). The trace for V = 22 
mV contains two-level fluctuations, as indicated by the two 
horizontal dashed lines. 



Fig. 1(a)]. The thickness of each layer is indicated in (. . .) 
in nanometers. The multilayer stacks are patterned into 
elliptic pillars with 160 x 60 nm dimensions by milling up 
to the middle of the PtMn layer. To crystallize CoFeB 
layers, the stack is annealed in 10, 000 Oe for 120 min at 
330 °C^2. 

All of the results presented here were obtained at low 
temperature (3-5 K) in the variable temperature insert 
(Oxford VTI) in a magnetic field (H) between -40, 000 
and 500 Oe. The schematic setup of our measurement 
system is shown in Fig. 1(a). A dc bias is applied under 
a constant current condition (I) by using the dc voltage 
source through a 100 kSl resistor. The differential resis- 
tance (dV/dl) and the dc bias voltage (V) are measured 
by a lock-in amplifier and a digital multimeter, respec- 
tively. Figure 1(b) shows the typical MR curve around 
the hysteresis loop of the free layer at 5 K. The clear 
square shape of the curve without any steps indicates 
that there is no pinning site on a macroscopic domain 
wall. The MR ratio defined as (R AP -Rp)/R P is 208%, 
where Rp and Rap are the sample resistances in the 
parallel (P) and antiparallel (AP) configurations, respec- 
tively. Figure 1(c) shows the typical bias dependence 
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of the differential resistance, where the solid and dashed 
curves correspond to the AP (—500 Oe) and P (500 Oc) 
configurations, respectively. In the differential resistance, 
as temperature decreases below 10 K, a peak structure 
appears around the zero bias in both P and AP config- 
urations, which is consistent with the several previous 
reportSi2Ai& Interestingly, in the P configuration, addi- 
tional satellite peaks seem to appear around V = —5 and 
5 mV. Such a feature may be related to the observation 
reported before^ 

To obtain the voltage noise spectral density Sy, we 
measure the time-domain voltage fluctuation signal by 
using a two-channel digitizer (National Instruments PCI- 
5922), which yields Sv by a fast Fourier transformation. 
In this process, two sets of voltage signals are simultane- 
ously measured after being independently amplified by 
two room-temperature amplifiers (NF Corporation LI- 
75 A). The cross-correlation technique is used here to re- 
duce the external noise and amplifier noise by long-time 
averaging. 52 The frequency range of our system is 100 
Hz to 200 kHz. In addition to Sv, the real-time voltage 
signal is also recorded by the digitizer at a sampling rate 
of 1 MHz. 

The Sv measurement is carefully calibrated with 
the thermal noise of several commercial resistors 
(MCY100R00T, MCY250R00T, MCY350R00T, and 
MCY1K0000T) with a precision of 0.01%. The typi- 
cal resolution of Sv for shot noise estimation is below 
10~ 20 V 2 /Hz. As a result, we achieved an experimental 
precision for the Fano factor well below 1%. We mea- 
sured three devices with the same geometry (samples 1, 
2, and 3) which are made out of the single wafer, and 
obtained consistent results. 



B. Analysis of noise 

In the present experiment, we found that the voltage 
noises in MTJs consist of a white noise (S'whito), which 
is a frequency-independent component, 1// noise (Si//), 
and RTN (Srtn)^ 



Sv — S w hitc + Si 



If 
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The white noise (S w hite) is attributed to the thermal agi- 
tation of electrons (thermal noise) and the partition pro- 
cess of electrons (shot noise). Then S'whito is described 
by 
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where k B is Boltzmann's constant, e is the electron 
charge, Rd is the differential resistance (dV/dl) at a given 
V (or I), and F is the Fano factor. The 1// noise (Si//) 
in MTJs is parameterized by 2 ^ 



Af ' 



where a*, f, and A are the Hooge parameter, fre- 
quency, and junction area, respectively. The RTN is the 
fluctuation between two levels, where Srtn exhibits a 
Lorcntzian character in an ideal case. 

Typical results of Sv for the AP configurations for 
V = 12, 22, and 55 mV at 5 K are shown in Fig. 2(a). In 
the data, the resistor-capacitor (RC) damping owing to 
the capacitance (760 pF) of the measurement lines has 
already been corrected for^ In the inset of Fig. 2(a), 
which shows the region between 160 and 200 kHz, the 
spectra are almost flat, and thus the increase of Sv with 
increasing V corresponds to the shot noise. In contrast, 
the spectra for the low-frequency region strongly depend 
on the frequency, as shown in the main panel of Fig. 2(a). 
The spectra for V — 55 and 12 mV are clearly dominated 
by the 1// noise except for the high-frequency region. 
The spectrum for V — 22 mV has a clear Lorentzian char- 
acter^ indicating that the source of the noise is two- level 
fluctuation of the resistance, namely, RTN. Figure 2(b) 
shows the real-time voltage signal for V — 12, 22, and 
55 mV. Apparently, the time-domain signal has a two- 
level nature only for V = 22 mV, as indicated by the two 
horizontal dashed lines in Fig. 2(b). 

We performed a histogram analysis between 160 and 
200 kHz to estimate Swhite^ To investigate the 1 // noise 
and RTN, we first subtract the white noise component 
from the measured Sv and then define the Hooge pa- 
rameter a as 
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where the frequencies fi and fi are chosen to be 1 and 
100 kHz, respectively. The value of a thus obtained 
equals the a* in Eq. (3) when RTN is negligibly small. 
In this study, we mainly use the a that is a well-defined 
parameter even if RTN is present. 



III. RESULTS AND DISCUSSION 

A. Frequency-dependent noise 

We start with the experimental result of the frequency- 
dependent noise. There have been several studies on the 
field dependence of the 1// noise in Al 2 C^-based 2 ^ 32 ! 34 
and MgO-based^l^^^i 5 -! 4 ! MTJs. In these reports, 
the measured 1// noise consists of field-independent and 
field-dependent components. The former component has 
a nonmagnetic origin (charge trap in the tunneling bar- 
rier), whereas the latter component has a magnetic ori- 
gin (magnetic fluctuations and domain wall motion in 
the free and/or fixed magnetic layers). Recent studies on 
the 1// noise in MgO-based MTJs have shown that the 
nonmagnetic 1// noise is no longer important owing to 
the improvement of tunneling barrier quality. 43 For RTN, 
both magnetic and nonmagnetic origins were reported 
in AlsOs-based 2 ^ 3 ^ and MgO-based 3 ^^^ MTJs. 
Nevertheless, there are only a few systematic studies on 
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FIG. 3. The magnetic field and bias current dependence of a 
measured at 5 K for sample 1. (a) MR curves of this MTJ and 
\dR/dH\/ R values for P (dashed line) and AP (solid line) con- 
figurations numerically estimated from the MR curves. The 
vertical dashed line marks H = H c at which \dR/dH\/R for 
AP and P have the same value, (b) The measured a for AP 
(solid mark) configuration, where we marked the peaks of a 
as A and C (upper panel) , and a color plot of a as a function 
of magnetic field and bias current (lower panel). The dashed 
line corresponds to a in the upper panel, (c) Counterpart of 
Fig. 3(b) for the P configuration, where the peaks of a are 
marked as B and D. a values for the two configurations were 
measured in almost the same bias voltage region. 



the magnetic field and bias voltage dependence of RTN. 
Here we focus on the magnetic 1// noise and RTN in the 
hysteresis loop of the free layer. 

Figure 3(a) shows the resistance and \dR/dH\/R at 5 
K, which is obtained as the field is ramped from —500 to 
200 Oe for AP and from 500 to -240 Oe for P configura- 
tions. The value of a measured at each field for AP and P 
configurations is shown in the upper panels of Figs. 3(b) 
and 3(c), respectively. It is clear that the magnetic field 
dependencies between a and \dR/dH\/R resemble each 
other for the two configurations: 54 They increase with 
the reversal of the free layer, except for a few peaks ob- 
served in a [peaks A, B, C, and D in Figs. 3(b) and 3(c)]. 
Whereas the measured spectra are dominated by the 1 // 
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FIG. 4. Field sweep rate dependence of dc resistance. The 
solid line is the dc resistance measured with the slowest sweep 
rate (50 Oe/min) as the field was ramped from —500 to 240 
Oe for AP and 500 to —300 Oe for P. Results for other sweep 
rates (200, 500, and 1000 Oe/min) are plotted as differences 
from the slowest one. 



noise outside of these peak regions, strong enhancement 
of Lorentzian character is always seen in such peak re- 
gions. The lower panels of Figs. 3(b) and 3(c) show image 
plots of a as a function of the magnetic field and the bias 
current for the AP and P configurations, respectively. In 
these figures, the ranges of the bias current for the two 
configurations are set to be almost the same with respect 
to the bias voltage. The curves shown in the upper panel 
of Figs. 3(b) and 3(c) correspond to the cross sections 
indicated by the dashed line in each image plot. 

We define H c as the field at which \dR/dH\/R gives 
the same value for the AP and P configurations [see 
Fig. 3(a)]. The behaviors of a and \dR/dH\/R are found 
to be well symmetric with respect to H c , except for the 
peaks. The values of a for AP and P configurations at 
their baselines are 3.3 x 10 -12 and 2.7 x 10~ 13 /mi 2 , re- 
spectively. This subtle difference is presumably due to 
the magnetic fluctuation of the fixed layer. Remarkably, 
the peak positions are roughly symmetric with respect 
to H c (namely, peaks A and B, and peaks C and D). Fi- 
nally, these symmetric behaviors with respect to H c are 
always observed on all three samples. We will discuss 
the implications of these observations later in Sec. HID 
again. 

The influence of the magnetic field sweep rate on the 
MR curve is shown in Fig. 4, where H c is also a relevant 
parameter. The resistance measured by each sweep rate 
is plotted as the difference of the data from that obtained 
with the slowest sweep rate (50 Oe/min), where the field 
was ramped from —500 to 240 Oe for AP and from 500 
to —300 Oe for P configurations. Although initially there 
is no sweep rate dependence on the resistance for either 
configuration, a sweep rate dependence appears as the 
magnetic field crosses H c . This observation strongly in- 
dicates that the free electrode feels zero effective field at 
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FIG. 5. a at several bias currents vs \dR/dH\/R for P (open 
symbols) and AP (solid symbols). The inset shows the full 
scale of the same data for AP. The dashed line is a guide to 
the eye. 



H c . The shift of H c from zero field is explained by taking 
a magnetic interaction between the free and fixed layers 
into account. 



where S m , /j,q, and x'm denote, respectively, the spectral 
density of magnetic fluctuation, the vacuum permeabil- 
ity, and the imaginary part of the magnetic susceptibility. 
Then by using a typical equation of 1// noise [Eq. (3)] 
and the Kramers-Kronig relation, Eq. (5) reduces to 3 ^ 
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where 2m and Ai? are the respective changes of the mag- 
netic moment and resistance associated with the reversal 
of the free layer. Thus, the measured linear relation be- 
tween a and \dR/dH\/R for the P state in Fig. 5 can 
be explained by Eq. (6); namely, the origin of 1// noise 
is thermal agitation of the magnetic moment of the free 
electrode. The dashed line corresponds to Eq. (6), where 
we take the saturation magnetization Ms of CoFeB to 
be 1.3 x 10 6 A/m.— From rough estimation, we obtained 
fmax/fmin ~ 10 10±2 , which is consistent with the result 
fmax/fmin ~ 10 9 that Ingvarsson et al£2- reported. The 
nonlinear behavior with \dR/dH\/R beyond ~1 mOe -1 
is possibly caused by deviation from the thermal equi- 
librium state of the free layer. In fact, the measured 
dc resistance corresponding to this region strongly de- 
pends on the field sweep rate (see Fig. 4). According to 
previous workS ) 43 ' 55 i 56 the linear relation between a and 
\dR/dH\/R corresponds to constant magnetic losses and 
has been shown by Stearrett et al£^ in the P and AP 
configurations. 



B. 1// noise 

Now, we discuss what we can learn from the 1/f 
noise. Previously, Ingvarsson et al£Q scaled a ver- 
sus \dR/dH\/R and explained the result in terms of 
a fluctuation-dissipation (FD) relation. The scaling 
of the 1/f noise by the FD relation has been tested 
against Al 2 3 -based22i2I and MgO-based^i MTJs with 
micron-sized junction areas. Here, we test the scaling by 
the FD relation on the free layer of MgO-based MTJs 
with submicron-sized junctions following Ingvarsson et al. 
We replot the data in Figs. 3(c) and 3(d) and obtain the 
relation between a and \dR/dH\/R as shown in Fig. 5. 
Remarkably, the values of a for the P configuration (open 
symbols) are almost proportional to \dR/dH\/R. For the 
AP configuration, the a values (solid symbols) first de- 
viate from the dashed line (drawn to guide the eye) be- 
cause of peaks A and C, and then fall onto the dashed 
line again. In the inset of Fig. 5, a for the AP configura- 
tion is shown over the wider range of \dR/dH\/R. After 
dropping onto the dashed line, a for AP rapidly increases 
with \dR/dH\/R beyond -1 mOe -1 . In this region, the 
measured spectra exhibit a complicated enhancement of 
Lorentzian character possibly from several RTNs. 

If one assumes thermal equilibrium, then the FD rela- 
tion for this magnetic system is given by 30 

Sm = ^xl(f), (5) 



C. Random telegraph noise 

Previously magnetic RTN have been observed in 
AlaOa-based 3 ^ and MgO-based 3 ^ 3 ^ MTJs, which 
are typically sensitive to the magnetic field and bias volt- 
age. Although it is difficult to deal with the RTN system- 
atically in general, a few authors successfully estimated 
the effective magnetic moment of the fluctuator and dis- 
cussed possible origins of magnetic RTN. In AlaO^-based 
MTJs with micron-sized junctions, Ingvarsson et al 3 ® 
and Jiang et al. 3 - suggested a small rotation of a single 
domain or domain wall hopping between pinning sites. 
Recently, Herranz et al.~ reported that the magnetic 
RTN is caused by magnetic inhomogeneities and domain 
walls in free and fixed layers. 

We observed a strong enhancement of the Lorentzian 
components in regions with specific bias voltages and 
magnetic fields. Here we consider peaks A, B, C, and D in 
Fig. 3(b) as typical examples. Generally, the Lorentzian 
component in the noise spectral density indicates that 
the noise is caused by the fluctuation between two levels, 
namely, RTN. The real-time voltage signals near peak A 
for H = —300, —320, and —360 Oe are shown in the left 
panel of Fig. 6(a); the right panel shows the correspond- 
ing histograms for the voltage signals (2 x 10 5 points for 
200 ms). These signals have two distinct voltage levels, 
and hence we attribute the measured Lorentzian compo- 
nents to RTN. The histograms are fitted by a double- 
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FIG. 6. (a) Magnetic field dependence of the voltage fluctua- 
tions (offset for clarity) measured around peak A in Fig. 3(b) 
(left panel) and the corresponding histograms (dots) (2 x 10 
points) (right panel). The solid and dashed curves are the fit- 
ted curves (double Gauss function) and each Gaussian com- 
ponent, (b) Logarithm of the rate of the dwell times for 
I = —42 fiA (open circle) and 42 (j,A (solid circle) as a func- 
tion of H. The dashed and solid lines are the results of the 
linear fitting, (c) Bias dependence of estimated fluctuation 
magnetic moments. 



where E, Am, and 1/tq are the field-independent acti- 
vation energy, the total magnetic moment of the fluctu- 
ator, and an attempt frequency, respectively (where we 
take + for the 1 state and — for the 2 state). To esti- 
mate the effective magnetic moment parallel to H (Amy), 
we take the ratio of each dwell time : 30 ' 31 ' 46 t\/t2 oc 
exp (— 2Am||iJ/fcsT) . By a linear fit ofln(ri/r2) versus 
H [see Fig. 6(b)], we estimated the effective magnetic 
moment Amy to be 1.6 x 10 3 /is for both positive and 
negative currents (I = ±42 fiA), where hb is the Bohr 
magneton. This value corresponds to 0.08% of the total 
magnetic moment of the free electrode (2.0 x 10 6 /is). 
However, the change in resistance Ai? resulting from 
RTN was estimated from the real-time voltage signal to 
be 0.13 f2, which is about 0.035% of the total resistance 
change for the reversal of the free layer, which is con- 
sistent with the above value (0.08%). This observation 
strongly suggests that the observed RTN is responsible 
for the magnetic fluctuator. The bias dependence of the 
estimated Amy is shown in Fig. 6(c). With increasing 
bias current, Amy is monotonically decreased for both 
bias polarities. This means that as the injection current 
gets higher, smaller magnetic fluctuations are allowed. 

The Amy values in the previous report s 30 ' 31 ' 46 are 
larger than those in our result by at least two orders 
of magnitude. Here we discuss the origin of the fluctu- 
ator that contributes to the measured magnetic RTN. If 
one assumes full reversal of a single domain in the free 
layer, a typical area size of the magnetic fluctuator is 
5.9 x 10 -6 /jm 2 , where the junction area is 7.4 x 10~ 3 /jm 2 . 
This small value and the absence of any step in the mag- 
netic hysteresis loops in Fig. 4 indicate that the contri- 
bution of a macroscopic domain wall can be ruled out. 
The estimated Amy and Ai? for the AP configuration 
have almost the same percentages against the full rever- 
sal of the free layer. Finally, enhancement of magnetic 
RTN is observed for both configurations. Based on these 
results, we attribute the fluctuator to two quasistablc 
single-domain states with some strain in the free layer 3 s - 



Gaussian function, which is shown in the same panel by 
the dashed lines. Let us call the two states "1" and "2." 
From the area of each Gaussian, which is proportional to 
the dwell time of each state (n and T2), we estimate the 
ratio of the dwell times between these states (T1/V2). Fig- 
ure 6(b) shows a logarithmic plot of T\jr% versus H for 
/ = ±42 fiA. Remarkably, logTi/r 2 is proportional to the 
magnetic field. The strong dependence of the dwell time 
for each state as a function of magnetic field indicates 
that the RTN is due to a magnetic fluctuator. 

We assume that two states with dwell times of t\ and 
T2 have activation energies of E ± Am • H. By further 
assuming the Arrhenius relation, we express the dwell 
times Ti and Ti as 

1 1 / £±Am-H\ 

^)" eXP l k^)' (7) 



D. Crossover from the 1/f noise to RTN 

To understand the excitation mechanism of the RTN 
we focus on the bias dependence of the noise property 
at peak A. Figure 7 represents the frequency-dependent 
component of measured spectra at —300 Oe (AP config- 
uration) for V = 12, 23, 30, and 50 mV. Although pure 
1/f noise is observed in the low-bias region as shown in 
Fig. 7(a), RTN is dominant at V = 23 mV, resulting 
in the strong enhancement of the Lorentzian component 
("RTN1") [Fig. 7(b)]. The characteristic frequency of 
the RTN is given by the full width at half maximum of 
the Lorentzian (/ c ). For example, f c = 20 kHz for RTN1 
in Fig. 7(b). With increasing bias voltage, the character- 
istic frequency of RTN 1 increases and another Lorentzian 
component ("RTN2") with different f c appears [Fig. 
7(c)]. Finally, the 1/f noise becomes dominant again 
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FIG. 7. (a)-(d) The frequency-dependent component of the 
spectrum for V = 12, 23, 30, and 50 mV, respectively. The 
solid curves are the fitted curves obtained by summing the 1/f 
noise, RTN1, and RTN2, where each component is shown as 
a dashed curve. 

in the spectrum [Fig. 7(d)]. To fit the obtained spectra 
as a summation of 1/f noise (Sim) and a few Lorentzians 

(SrtnO, we use PIP and Ptl\ l + (///ci) 2 ], respectively. 
Here 7 and fd are the spectral exponent and character- 
istic frequency of the Lorentzian for "RTNi" (i = 1 and 
2). Outside the peak region, the estimated 7 is close to 
1, whereas it is reduced from 1 as the Lorentzian compo- 
nents are enhanced. By taking the Dutta-Dimon-Horn 
models in which the 1/f noise is the result of a su- 
perposition of many RTNs with a broad distribution of 
activation energies, this suppression of 7 accompanied by 
the enhancement of RTN can be explained by a change 
of the energy distribution of the magnetic fluctuators. 

The bias dependence of the parameters characterizing 
RTN at —300 Oe for the AP configuration is summarized 
in Fig. 8. V c ~ 21 mV is the bias voltage where a 
shows its maximum. The effective magnetic moment of 
the fluctuator, Amy, and the change in resistance, AR, 
are compared in Fig. 8(b). The behaviors of a, Am.11, 
and AR are similar to each other in that they have their 
maxima at V c and then decrease as the bias is increased. 
Finally, both Amy and AR become almost half of their 
maxima at the bias at which the peak of a disappears. 



Figure 8(c) shows that, although the estimated f c for 
RTN1 and RTN2 exhibits almost an exponential increase 
with the bias for their whole curves, there is a kneelike 
structure in the curve for RTN1 at V c . We fit f c for RTN1 
to f c = fo exp (— E/ksT*) using the Arrhenius relation 
and taking Joule heating into account, where fo is an 
attempt frequency. We also assume that the activation 
energy E is independent of the bias voltage and that the 
effective temperature T* = T + Cy 2 . The solid curve 
in Fig. 8(c) corresponds to the fitted curve, where the 
estimated E and ( are 1.3 x 10" 21 J and 1.3 x 10 4 K/V 2 , 
respectively. Comparing this result to the previous one^ 
in which the activation energy for nonmagnetic RTN in 
Al 2 03-based MTJs was found to be 0.3 x 10 -19 J, we see 
that our value is much smaller. This is consistent with 
the absence of nonmagnetic RTN in high-quality MTJs 
in the low-temperature, low-bias regime. 

Returning to the magnetic configuration dependence 
of the RTN, one sees that a values for peaks A and B are 
enhanced at roughly symmetric positions with respect to 
H c [see lower panels in Figs. 3(b) and 3(c)]. The values 
of a for peaks A and B have their maxima at almost the 
same bias voltage, whereas the amplitude of the maxi- 
mum for peak A is almost seven times as large as that 
of peak B. We assume that the same magnetic fluctuator 
contributes to these two peaks. The symmetric hysteresis 
loop of the free layer with respect to H c also supports this 
assumption (see Fig. 4). Here, the difference in the peak 
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amplitude between A and B cannot be accounted for by 
conventional Joule heating. That is, because the Joule 
heating effect is proportional to the bias current at a fixed 
bias voltage, this effect for peak A is only 1/3 that for 
peak B. These results may imply a spin-dependent heat- 
ing process of a localized spin system in the free layer 
depending on the magnetic configurations. 



E. Shot noise 

Finally, we discuss the shot noise to connect our pre- 
vious work with the present one and to support fur- 
ther evidence of our claim made before. The frequency- 
independent component of the spectrum is well described 
by Eq. (1). F describes how the noise deviates from 
the Poissonian value and thus characterizes the partition 
process of the electron tunneling. Conventional tunnel 
junctions exhibit F — 1, reflecting the Poissonian pro- 
cess. Regarding Al 2 03-based MTJs, after the first report 
of the full shot noise (F ~ l)r^ reduced Fano factors 
ranging from 0.45 to 1 were reported] 15 ! 16 This reduction 
can be explained by the sequential tunneling modelji 6 - 
where the process of two-step tunneling through impu- 



rities within the barrier is assumed. In this theoretical 
model, F strongly depends on the asymmetry of the each 
tunneling and can be 0.5 to 1. 

Full shot noise in MgO-based MTJs was reported by 
Guerrero et al. and some of the authors of the present 
paper in experiments where the MgO barriers were as 
thick as 3 and 1.5 nm, respectively 1 7 ' 18 This indicates 
that MgO-based MTJs are free from the process through 
impurity sites, possibly resulting from the high quality 
of the crystallized MgO. Later, we reported the suppres- 
sion of the Fano factor for the P configuration (typi- 
cally 0.91) with a 1.05-nm-thick MgO barrier, whereas 
the F for the AP configuration is almost unitary (typi- 
cally 0.99) To explain this subtle reduction of F within 
the above sequential tunneling model, we have to assume 
a very asymmetric barrier (1:100), which is unrealistic as 
the barrier of our MTJ is thin. Moreover, by using this 
model, we cannot explain the magnetic configuration de- 
pendence of the reduced F, and hence we can rule out 
this scenario. We note that all of the experimental re- 
sults for F values including the thick-barrier case can 
be explained by recent theoretical work based on a first- 
principles calculation^ Namely, the Fano factor, which 
is reduced from 1, is a direct consequence of the Pauli ex- 
clusion principle, signaling that there are coherent chan- 
nels with high transmission probabilities through the epi- 
taxial MgO barriers. Thus, our result on the shot noise 
gives unique evidence for coherent tunneling through a 
crystallized MgO barrier. 

We show in Fig. 9 the typical results of the white-noise 
component Sy and a at 3 K as a function of bias voltage. 
In Fig. 9(a), S v is nicely fitted by Eq. (2) for \V\ < 10 
mV. For \V\ > 10 mV, a increases due to enhancement 
of RTN for both configurations, and, as a result, Sv de- 
viates from the fitted curve. It is noted that the effect of 
the 1/ / noise on Sv is negligibly small in this bias range. 
The above observation indicates that to accurately esti- 
mate the Fano factor immune to the frequency-dependent 
noise, it is necessary to analyze the noise in the low-bias 
regime at low temperature. Experimentally, we can sup- 
press magnetic RTN by stabilizing the magnetization by 
applying a large magnetic field for the P configuration. 
In fact, in Fig. 9(b), the enhancement of RTN is sup- 
pressed, and thus Sy is almost perfectly fitted by the 
shot-noise formula with F — 0.92. 



IV. CONCLUSION 

The low-frequency noise properties in submicron-sized 
CoFeB/MgO/CoFeB-based MTJs with thin tunneling 
barriers are systematically investigated. The measure- 
ments are carried out at low temperature by controlling 
the magnetic field and bias current (voltage), where we 
focus on the noise properties on the magnetic hystere- 
sis loops of the free layer. A clear correlation between 
a and \dR/dH\/R is observed. The nice scaling of the 
observation by the FD relation indicates that the main 
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origin for the 1// noise is thermal magnetic fluctuation 
of the free layer. RTN is observed inside the magnetic 
hysteresis loops for both configurations. We found that 
this is due to the magnetic fluctuation between two qua- 
sistable single-domain states with some strain in the free 
layer. Although the 1/f noise is almost independent of 
the magnetic configuration, RTN is remarkably enhanced 
for the AP configuration and at specific bias voltages. 
Such results indicate a spin-dependent activation process 
of RTN. Shot noise measurement gives us quantitative in- 
formation for coherent tunneling. 

Our study shows that a systematic study on the noise 



of MT Js is possible by using devices with well-crystallized 
thin MgO barriers and by measuring the noise at low tem- 
perature and low bias. Further study of the noise proper- 
ties in MgO-based MTJs with various barrier thicknesses 
is necessary to systematically elucidate the mechanism of 
the noise and improve the device properties. 
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